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Abstract 

We present a one dimensional model for diffusion on a hierarchical tree structure. 
It is shown that this model exhibits aging phenomena although no disorder is present. 
The origin of aging in this model is therefore the hierarchical structure of phase space. 

1 Introduction 

Strongly disordered systems have been a focal point of research in recent years. Among 
the most studied materials are magnetic systems with impurities and especially spin glasses 
0. Their dynamics can be characterized by glassy features which are also encountered 
in many different systems such as structural glasses Q, polymers in random potential ||, 
protein folding [|J, dirty superconductors 0, charge density waves with impurities and 
also areas like biological evolution 0. 

Spin glasses are magnetic materials with structural disorder, e. g. alloys of a magnetic 
and a non-magnetic metal as for example AgMn and CuMn. An overview of different spin 
glass materials can be found in 0. In the beginning, twenty years ago, the main interest 
was to find the equilibrium properties of spin glasses, since it was expected that a physical 
system reaches it's equilibrium in finite time. In spin glasses, however, the dynamics is 
governed by non-equilibrium properties such as aging which was found almost ten years 
later [|]. 

In this context aging describes the striking effect that magnetic properties in spin 
glasses depend drastically on their history (or age) in the frozen phase. Although it is 
mainly referred to aging in the context of spin glasses it can be observed in other disordered 
substances || also. 

Up to now it is still a major task in theoretical physics to understand the origin of 
aging, although many different models have been proposed in recent years. Most of these 
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are phenomenological theories like the droplet model |I(J, the domain growth theory [JTT 



or the trap model A different approach is given by mean- field models of Ising spin 
glasses [13] and models with a hierarchical structure in phase space [14]. Most of these 
models are rather successfull in fitting experimental data although some of them differ in 
their final conclusions concerning the scaling of time dependent properties. 

In order to discriminate between these differing scaling assumptions many experiments 
have been carried out ||15|| , but the situation remains unclear. Therefore numerical studies 
play an important role as one can obtain detailed information about the dynamical pro- 
cesses and the spatial correlations in glassy systems. Much progress has been obtained 
recently in the numerical studies on aging in spin glasses and other disordered systems 
In numerical investigations the autocorrelation function can be calculated directly 
and one finds a crossover from a slow quasi-equilibrium decay for t <C t w to a faster non- 
equilibrium decay for t 3> t w . The functional form of these decays, e. g. for the three 
dimensional Edwards- Anderson spin glass model, is algebraic and has the scaling form 



(1) 



with &t(u) = c y , for y = 0, and $r(y) oc y x ( T, )- x( - T \ for y — > oo JT^]. As an important result 
one gets two different exponents which are characteristic for the two dynamic regions. 
Since the algebraic decay of correlation functions is also known from diffusion models 



in ultrametric hierarchies |I^,|T8[, which can be solved analytically, and the analogy to 
the structure of Parisi's symmetry breaking scheme in the solution of the EA-model a 
hierarchical ansatz to explain aging phenomena seems rather natural. In recent years 
different hierarchical models for aging were proposed [14|. Most of these models incorporate 
randomly chosen distributions of diffusion rates and, as a consequence, random energy 
barriers and trapping times. With further assumptions concerning the functional form 
of these distributions one gets aging phenomena like the crossover from quasi- to non- 
equilibrium dynamics with an algebraic scaling form according to ([!]). 



Taking long range diffusion in an ultrametric space []I7| as starting point one can con- 
struct a hierarchical model which shows aging phenomena as a direct consequence of the 
structure and not of random properties. This model will be presented in the next section 
followed by numerical results and some conclusions. 



2 The Model 

In contrast to the known hierarchical models for aging the presented model incorporates 
no kind of disorder or randomly distributed energies. The main idea is to take into account 
the different types of dynamics found in spin glasses. At short time scales t <C t w one finds 
quasi-equilibrium dynamics and at long times t ^> t w non-equilibrium dynamics with a 
crossover region for times t ~ t w . These different types of dynamics are commonly related 
to regions in phase space which are accessible at different time scales. 
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Figure 1: A hierarchical tree structure with three types of states characterized by different 
line styles: thick solid lines represent the surrounding non-equilibrium tree, dashed lines 
represent the quasi-equilibrium subtrees and thin solid lines represent the attractor region. 
The corresponding branching numbers are n = 3, m = 2 and n a = 4; the height of the tree 
is k = 3. 



A rather natural ansatz to include these different types of states is to extend the diffu- 
sion on a hierarchical tree structure in which every state is equal by a more sophisticated 
tree structure which has at least two different types of states. The first type charactarizes 
the quasi equilibrium dynamics. These states are grouped as leaves in regularly branching 
trees of varying heights where standard long range diffusion in an ultrametric space takes 
place [fTTH - The second class of states is found in a regular tree structure in such a way that 
the different equilibrium trees are embedded in this surrounding tree, i. e. every branching 
point of the non-equilibrium tree is the root of a quasi-equilibrium subtree. The resulting 
structure is shown in Figure |l|, where a third class of states is included in order to represent 
some kind of attractor states for the equilibrium trees. The over-all height of the tree is 
k, which equals the number of different hierarchies in the tree structure. Figure [l] shows 
a tree of height k = 3 and the different trees can be recognized by the line styles, i. e. the 
thick solid lines represent the surrounding non-equilibrium tree with the branching number 
denoted by n (=3), dashed lines represent the quasi-equilibrium tree with the branching 
number m (= 2) and thin solid lines represent the attractor states with the branching 
number n a (= 4). 

The resulting diffusion rates in the master equation, defining the dynamics in the whole 
tree structure, are chosen to be unsymmetric which takes into account the different classes 
represented by the states. Let the system start at an arbitrary non-equilibrium state Qq. 
We then denote the probability to be found in the whole subtree of height I (1 < I < k) 
including Q with Qi, the probability that the system can be found in a state of an quasi- 
equilibrium subtree which can be reached from the initial state by crossing I hierarchies 
with Pi and finally the probability to be found in the attractor states corresponding to 
the equilibrium subtree Pi with Oi. The allowed long-range hops in this tree structure are 
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formulated in the following way ||19|| : 

1. From a state of the non-equilibrium type every other state in the tree can be reached 
with an rate with i being the number of hierarchies between the states. 

2. From an equilibrium state in Pi every state within Pi can be reached with a rate r; 
(1 < i < I) and outside of Pi every state in the tree outside from Qi is allowed with 
the same diffusion rate as for the non-equilibrium states. 

3. From the attractor states Oi hops to every state in the tree are allowed, but with the 
rate si into the corresponding equilibrium tree Pi and with the rate T{ {I < i < k) 
into all other states. 

These rules include that the system has to cross one hierarchy more in order to leave 
an equilibrium subtree because the other states within Qi are excluded by the rules. This 
asymmetry is chosen to generate the equilibrium subtrees Pi as dynamical traps with 
characteristic escape times corresponding to their height. Using these dynamical rules the 
master equations can be formulated in a closed way 

^-Pi{t) = -AiP^ + BiQi^ + BiiQi^-P^-Oi^n^ + BiSiOiit), 
at 

^-OW = -AiOi^ + ani^-BiSiOi^-iAi-Bi-ayi^Oiit) (2) 
at 

+c l (Q l (t)-P l {t)-o l (t))r l „ 1 , 

^-Qi(t) = -At Qi(t) + AiQi(t) . 
at 

To write the equations in a compact form the following abbreviations were used: 

k-l 

At = ^(A m -^)r,, (3) 
i=i 

k-l 

Vi(t) = ^(Q <+1 (t)-Q < (t)-P <+1 (t))r < . 

i=i 

The constants Ai, Bi and C; denote the number of states in a whole subtree Qi, an equi- 
librium subtree Pi and the attractor states Oi respectively. The diffusion rates ri and si 
are as mentioned above. 

The master equations can be represented by a matrix which is an upper triangular 
matrix despite some entries on the lower next- diagonal. In principle this system can be 
solved for arbitrary numbers of states and diffusion rates. In order to refer to some exact 



results known for the "classical" diffusion models |T7[ the branching numbers of the states 



are taken to be constant. Therefore the number of states as defined above are: 



n" — m" n a l — n l 



Ai = n l + m h n a 



n — m n a — n ' 



Bi = m\ (4) 
Ci = n n : 
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with the branching numbers n, m and n a for the non-equilibrium, equilibrium and attractor 
states respectively. If the ratio of the diffusion rates is chosen to be constant one simply 
gets r\ = r l and si = s . With this choice the master equations are solved numerically and 
some of the results are presented in the following section. Analytic calculations concerning 
the exact solution of the model can be found elsewhere |T9|j . 



3 Results 

Aging phenomena are common for magnetization curves in real spin glass materials. 
The corresponding function describing the main dynamical effects for spin- or diffusion- 
models is the autocorrelation function. As for magnetization experiments it is also possible 
to perform waiting time dependent calculations numerically. For such purposes one has 
to define the correct waiting time dependent autocorrelation function C(t,t w ) which is in 
case of this hierarchical model given by 

k k k k k 

C(t, t w ) = Pjo(t) c, a Pi(t w , i) + Q (t) c i0 Q {t w , i)+y^ y O l0 (t) ^ c u Oi(t w , i) , 

1=1 i=l i=0 1=1 i=l 

(5) 

where Pi(t w ,i) is the probability that the system has reached an equilibrium subtree of 
height I with i hierarchies between this subtree and the starting site after the waiting time 
t w has elapsed. The constant cu denotes the number of such subtrees to be found. The 
function Pio(t) describes the probability that the system is still in the same state or has 
returned to it after the additional time t has elapsed as it was in at the end of the waiting 
time t w . The probability functions of the other states are denoted in the same way. 

The waiting time dependent autocorrelation function C(t,t w ) as defined in (|^) is the 
focal point of interest in the following. The system parameters are the branching numbers 
m, n, n a and the diffusion rates r and s. These parameters can be varyied in order to 
obtain numerical data. 

Figure [| shows the waiting time dependent autocorrelation function C(t, t w ) for param- 
eter values as given in the caption. A crossover from slow dynamics for small times to a 
faster decay at large times charactestic for aging can be seen very clearly. The crossover 
region is located at times t ~ t w . The decay obeys an algebraic time dependence although 



some oscillations characteristic for self-similar systems jl7j occur, especially in Figure 0(a). 
These oscillations show a strong effect on the decay exponent which can be calculated via 
the logarithmic time derivative of the autocorrelation function 

d\nC(t,t w ) 

7 = aEi— (6) 

The resulting exponent 7 of an algebraic decay C(t) oc t~~ f is shown in Figure [| The 
strong effect of the oscillations of the autocorrelation function on the exponent can be seen 
in part (a) of Figure |3|. As shown in part (b) the mean value 70 of the exponent taken 
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(b) 

,miiiii| Miiiiij iiinii| Miniii| iiniiii| iMiiini iMiiini 




Figure 2: The waiting time dependent autocorrelation function C(t,t w ) versus time t for 
varius t w . The system parameters are as follows: m = 3,n = 5,n a = 7, r = 0.02, s = 
0.3; t w = 10 4 ,10 7 , 10 10 ,10 13 in (a) and m = 2, n = 3, n a = 9,r = 0.04, s = 0.5; t w = 
10 2 ,10 4 ,10 6 ,10 8 in (b). 




Figure 3: The exponent 7 as defined in equation (|) for the autocorrelation function C(t, t w ) 
using the same parameter values as in Figure @(a). The left diagram shows the exponent 
with the strong oscillations and the right diagram shows the mean value 70 of the exponent 
calculated over one period of oscillation. 
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Figure 4: The exponents 7 e and 7 n versus the diffusion rate r for different branching num- 
bers m. The values of m are marked with different symbols: m = 2 (x),m = 3 (D),m = 
4 (o), m = 5 (*), m = 6 (o), m = 7 (A). The values of t w and s vary, n = 7 and n a = 9. 



over one period is constant for the two different time regimes t <C t w and t ^> t w . So the 
autocorrelation function satisfies the scaling relation (|IJ) with temperature independent 
exponents j e for t <^ t w and 7„ for t^> t w . 

Since the branching numbers and the diffusion rates play an important role in this 
model the exponents will at least depend on some of these parameters. The dynamics of 
this model is chosen to be close to the model solved in [|I7[] so it is expected that the values 
of the exponents 7 e and 7„ will be related to the decay exponent calculated analytically 
for the simple diffusion model [TH], i. e. 

InM 1 
r 7^|-lnM fOT R< Ji- < 7 > 

where M is the branching number in the tree and R the diffusion rate for crossing one 
hierarchy. The relation R M < 1 leads to an algebraic decay of the autocorrelation function, 
while one gets an exponential decay in the other case. 

The exponents j e and 7„ for different branching numbers m and diffusion rates r are 
shown in Figure |j. The dotted lines in diagram (a) mark the equilibrium exponent V 
according to equation (|7|). In the right diagram the 7-axis is shifted by one and the lines 
mark the values of T + 1. In both diagrams the numerical values are in excellent agreement 
with the theoretical prediction. As a consequence the values of n, n a and s are irrelevant 
for the dynamical behaviour of this model. This aspect and it's implications will be further 
discussed in the next section. 
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4 Conclusion 



A diffusion model on a hierarchical tree structure without any disorder was introduced. 
It was shown that the waiting time dependent autocorrelation function C(t,t w ) shows the 
characteristic dependence of the waiting time t w known as aging in spin glasses. Since the 
decay of the autocorrelation function is algebraic the scaling law ([!]) is satisfied with two 
constant exponents. The influence of the parameters on the two exponents were investi- 
gated and an excellent agreement with the equilibrium exponent F was found for short 
times t <C t w . After the crossover region a non-equilibrium exponent, in good agreement 
with r + 1, was found. The dynamic interpretation can be subsumed in the following way: 

• during the waiting time t w the system gets "trapped" in an equilibrium tree of height 
I corresponding to a trapping time r ~ t w \ 

• starting the measurement of C(t, t w ) after the waiting time causes the system to stay 
in this equilibrium subtree of height I; 

• the short time dynamics t <C t w is therefore characterized by equilibrium diffusion 
and the equilibrium exponent T is found; 

• in the crossover region t ~ t w the system slowly escapes out of this trap and the 
resulting exponent rises; 

• the long time behaviour t 3> t w corresponds to non-equilibrium dynamics with an 
exponent T + 1 ; 

• the attractor states 0\ together with a large diffusion rate s forces the system into 
the equilibrium subtrees; 

• since the system can be found in states of the equilibrium subtrees most of the times 
the dynamics is completely governed by the branching number of the equilibrium 
states. 

The presented model includes no kind of disorder but solely a hierarchical structure of 
phase space so it is obvious that the hierarchical organization of traps or valleys in the 
energy-landscape plays an essential role in explaining the aging phenomena. The results 
presented here are in good agreement with the scaling assumption ([!]) and coincide with 



the results of other hierarchical models for aging [13] and numerical investigations . 

Since the world of spin glasses and aging is much richer than discussed here, there are 
of course some additional phenomena which should be examined. In relation to experi- 
mental results one should consider temperature steps or more sophisticated temperature 
cycles. The model defined above is not directly temperature dependent, but with a straight 
forward relation between temperature and diffusion rates one can map temperature steps 
on corresponding steps in the diffusion rates [I5|. This should be the focus for further 
investigations. 
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